A&A manuscript no. 

(will be inserted by hand later ) 



Your thesaurus codes are: 

06 (08.16.5; 08.03.4; 02.12.3; 08.01.2; 08.13.1; 08.09.2 SU Aur) 



ASTRONOMY 

AND 
ASTROPHYSICS 



o: 
o. 

O' 
CN; 

Or 



> 

00 

o 
o> 
o 
o 
o 

Oh! 

6 

CO 

c3 



X 



Magnetospheric Accretion and Winds on the T Tauri Star 
SU Aurigae: 

Multi-Spectral Line Variability and Cross-Correlation Analysis/ 

J.M. Oliveira 1 2 , B.H. Foing 1 3 , J.Th. van Loon 4 , and Y.C. Unruh 5 

1 ESA Space Science Department, ESTEC/SCI-SO, P.O. Box 299, NL-2200 AG Noordwijk, The Netherlands 
e-mail: joana@so.estec.esa.nl, bfoing@estec.esa.nl 

2 Centro de Astrofisica da Universidade do Porto, Rua das Estrelas s/n, PT-4150 Porto, Portugal 

3 Institut d'Astrophysique Spatiale, CNRS/Univ. Paris XI, Bat. 121, Campus d'Orsay, F-91405 Orsay Cedex, France 

4 Institute of Astronomy, Madingley Road, Cambridge CB3 OHA, United Kingdom 
e-mail: jacco@ast.cam.ac.uk 

5 Institut fur Astronomie, Universitat Wien, Turkenschanzstr. 17, A- 1180 Vienna, Austria 
e-mail: ycu@astro.univie.ac.at 

Received 12 May 2000 / Accepted 21 August 2000 



Abstract. SU Aurigae is a T Tauri star that was well 
monitored during the MUSICOS 96 multi-site campaign. 
We present the results of the spectroscopic analysis of the 
circumstellar environment of this star, particularly of the 
Ha, H/3, Nai D and Hei D3 line profiles. The signatures 
of modulated outflows and mass accretion events are anal- 
ysed, as well as transient spectral features. We compute 
the cross-correlation function (CCF) of several pairs of 
(velocity bins in) spectral lines to better investigate the 
profiles' temporal variability. We found increasing time 
lags between the variability of Hei D3, Nai D and H/3. 
We propose this may be understood in terms of azimuthal 
distortion of the magnetic field lines due to the different 
rotation rates of the star and the disk. We find the slightly 
blueshifted absorption features in Ha, H/3 and Nai D to 
be anti-correlated with the accretion flow signatures. We 
propose that the transient absorption features in the blue 
wings of Ha, H/3 and Nai D (signatures of mass outflows), 
and flare brightenings are related to the disruption of dis- 
torted magnetospheric field lines. 

Key words: Stars: pre-main sequence - Circumstellar 
matter - Line: profiles - Stars: activity - Stars: magnetic 
fields - Stars: individual: SU Aur 
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* based on observations obtained during the MUSICOS 96 
campaign in which were involved: Isaac Newton Telescope 
(INT, La Palma), Observatoire de Haute Provence (OHP, 
France), McDonald Observatory (MDO, USA), Beijing Astro- 
nomical Observatory (BAO, Xinglong, China) and Canada- 
France-Hawaii Telescope (CFHT, Hawaii). 



1. Introduction 



1.1. Classical T Tauri Stars 



T Tauri stars (TTS) are low-mass (M < 3 M ) pre-main 
sequence stars which are normally younger than 3 x 10 6 



years (Appenzeller & Mundt 1989). They are of great in- 



terest as prototypes of young solar- type stars. Historically 
two sub-groups of TTS were defined, based on the strength 
of the Ha spectral line emission: classical T Tauri stars 
(CTTS) with strong Ha emission and weak (or "naked") 
T Tauri stars (WTTS). The enhanced activity of most 
CTTS as seen from the Balmer and metallic (permitted 
and forbidden) line emission is believed to be due to the 
presence of an accretion disk, whose UV and IR excess 
emission causes also broad energy distributions. 

Since the first studies of their line profiles it has been 
known that C TTS spectra show evidence of outflows (e.g. 
Herbig 1962). More recently also abundant observational 
evidence of funnel accretion flows onto CTTS has been 



published (e.g. Edwards et al. 1994). 

Konig l (|1991 ) proposed the application of the Ghosh 
& Lamb (1978) magnetospheric model, developed for ac- 
creting neutron stars, to explain the interaction between 
CTTS and their accretion disk. The accretion disk is dis- 
rupted at a few stellar radii from the central star (close to 
the co-rotation radius) by the presence of a stellar dipole 
magnetic field. According to Shu et al. (1994) and Paatz 
& Camenzind ( |1996 ), at the truncation point, the ionized 
disk material is loaded either onto inner closed magnetic 
field lines, accreting thereby onto the stellar surface, or 
it is loaded onto outer open magnetic field lines that can 
drive a disk- wind flow. 

Recently, some observational support for the magneto- 
spheric models was presented by Johns-Krull et al. ( 1999 ) 
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and Guenther et al. (1999) who measured kilogauss mag- 
netic fields at the surface of CTTS, enough to effectively 
disrupt the accretion disk at a few stellar radii. The as- 
sumption of a stellar dipole magnetic field aligned with 
the rotation axis completely determines the geometry and 
dynamics of the disk-star interaction. According to the 
magnetospheric models mentioned previously, the magnc- 
tosphere is axi-symmetric with respect to the rotational 
axis so no rotational modulation is expected in the wind 
and accretion signatures, unless the disk is inhomogeneous 
at the truncation radius. Furthermore, variations in the 
inflow and outflow signatures should be approximately in 
phase. But a stellar magnetic field that has solely a dipolar 
component is unlikely and this type of models has been 
criticised (Safier |199S| ; Alencar & Basri |2000f ). The evi- 
dence for magnetospheric accretion exists (hot spots and 
redshiftcd absorption features), but theoretical models are 
still far from explaining the variety of observations. 

1.2. SU Aurigae 

SU Aurigae has been classified as a CTTS (e.g. Bouvier 
et al. 1993 ), but also as the prototype of a separate class 
of T Tauri stars, the SU Aur type stars (Herbig & Bell 
198$ ). SU Aur has a G2 spectral type (Herbig |l960t ) , thus, 



due to the contrast with a hotter photosphere, only Ha 
and occasionally H/3, appear in emission in the spectrum. 
SU Aur has an exceptionally high projected rotational ve- 
locity of 60 km s- 1 (Johns-Krull |1996|; Unruh et al. |2000D 



and a rotational period of about 3 days. Photometric pe- 
riod measurements for SU Aur are very difficult to make. 



they observed in SU Aur. 

Using Ha, H/3, Nai D and Hei D3, we try to disentan- 
gle the two main contributions present in the line profiles, 
namely accretion and wind signatures. These lines probe 
different parts of the circumstellar environment. Thus, the 
temporal relation between the variability in the different 
lines depends on the geometry and physical conditions of 
the circumstellar environment, in particular of the inner 
disk and magnetosphere. For instance, variations in the 
accretion rate at the inner disk or outer magnetosphere 
will first affect low excitation energy features arising in 
this region. The higher energy lines that form in the inner 
accretion stream or in the accretion shock will be affected 
by this "perturbation" after a time lag related to the infall 
time through the magnetosphere. This delay is expe cted 
to be of the order of several hours (Smith et al. 1995 ). 

Previous SU Aur data sets typically had one to two 
spectra a night, even though over long time spans. Our 
data set, with finer time sampling, allows us for the first 
time to study the variations of the different components 
in spectral lines on short timescales. In the context of the 
magnetospheric model and using cross-correlation tech- 
niques, we have searched for time lags between the vari- 
ability of the mentioned spectral lines. In Sect. 3 we 
describe the implementation and results of the cross- 
correlation analysis of our data set. In Sect. 4 we present 
our interpretations in the frame of accretion and wind 
processes. In Sect. 5 we describe two transient absorption 
features observed in Ha, H/3 and Nai D and we analyse 
their velocity evolution. 



For instance, Bouvier et al. (1993) found multiple periods, 2. Description of the SU Aur Data Set 



with 2.78 days the most likely period. So far the best pe- 
riod estimates come from the analysis of the s pectr al vari- 
ations of the Balmer lines. Giampapa et al. ( 1993 ) found 
a period of 2.98 ± 0.4 days in the blue wing of the Ha 
line, later confirmed also in the blue and red wings of H/3 
by Johns & Basri ( 1995b ). Petrov et al. ( 1996 ) reported a 
period of 3.031 ± 0.003 days measured on the red wings 
of the Balmer lines. Our data set suggests a shorte r rota - 
tional period of about 2.6 to 2.8 days (Unruh et al. 2000] ). 



1.3. Diagnostics and Models of Magnetospheric Accretion 
From period analysis and equivalent width measurements 



of fitted components in Ha and H/3, Johns & Basri (1995b) 
observed that the H/3 redshifted absorption feature and 
the Ha blucshiftcd absorption feature (respectively ac- 
cretion and wind signatures) show a periodicity of about 
3 days and are approximately 180° out of phase in SU Aur. 
This led them to propose a generalization of the magneto- 
spheric model, in which the rotational axis and the mag- 
netic dipole axis are (slightly) misaligned: the misaligned 
"egg-beater" or oblique rotator model. It predicts that the 
signatures of mass accretion and disk winds should be ro- 
tationally modulated and approximately in anti-phase, as 



2.1. MUSIC OS 96 Campaign 

The MUSICOS (MUlti-SIte Continuous Spectroscopy) 
network was created to provide continuous spectroscopic 
coverage for stellar objects, whose variability timescales 
demand simultaneous coverage from several longitudes 
(Cata la & Foing [l990| ; Catala et al. |l993| ; Catala & Foing 
1994 ). To fulfill this aim, existing observational capabili- 
ties were used or dedicated instruments were developed. 
A MUSICOS campaign is organized every two years and 
several stellar targets are selected for each campaign. 

Our data set was obtained in November 1996 during 
the MUSICOS 96 multi-wavelength campaign, that in- 
volved five different observatories: Isaac Newton Telescope 
(INT, La Palma), Observatoire de Haute Provence (OHP, 
France), McDonald Observatory (MDO, USA), Bei- 
jing Astronomical Observatory (BAO, Xinglong, China) 
and Canada- France-Hawaii Telescope (CFHT, Hawaii). 
We obtained 126 echelle spectra of SU Aur, spanning 
10.5 days, with 3.5 days covered almost continuously. Ha, 
Nai D and Hei D3 were observed at all five observato- 
ries while H/3 could only be observed from the INT, OHP 
and CFHT. Thus, for this line, we obtained only 72 spec- 
tra. The data from the different observatories were re- 
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binned to a uniform resolution of 30 000, or approximately 
10 km s _1 . The journals of the observations and more de- 
tails about this campaign and data reduction can be found 
in Unruh ct al. (|2000). 



2. 2. Brief Summary of Variability in the Data Set 
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Fig. 1. The average (line plots) and variance (shaded 
plots) of (a) Ha, (b) H/3, (c) Nai D and (d) Hei D3. The 
y-axis scale on the left is for the average and on the right 
for the variance. For Ha and H/3 a rotationally broadened 
template profile is plotted (dot-dashed line in (a) and (b)) 
as a photospheric reference. The narrow feature that ap- 
pears in the variance profiles of Ha, H/3 and Nai D, be- 
tween — 150 and —50 km s _1 , is the signature of transient 
outflows (Sect. 5). The sharply peaked central feature vis- 
ible in the variance profile of Na I D is related to the nar- 
row interstellar component. The two components present 
in the Hei D3 profile are also plotted (dot-dashed line in 
(d)); the smaller features appearing at 180 and 230 km s _1 
are telluric water residues. The vertical dotted lines ap- 
pearing in all the graphs are the velocity bins chosen for 
the CCF analysis (Sect. 3.2); the binsize is 5 km s _1 . 

We highlight in this subsection the main characteris- 
tics of the variability in our data set. The period anal- 
ysis was performed by Unruh et al. ( 1998| 200C). They 
computed the periodograms of the above mentioned spec- 



tral lines (see figures in Unruh et al. 200C). We inde- 



pendently determined the following periods: 2.63 days 
in Hei D3 [—75:175] km s _1 and in the red wings of 
Nai D [30:150] km s _1 , and 2.80 days in the red wing 
of H/3 [50:175] km s _1 and in the far blue wing of Ha 



[—375: — 175] km s _1 . This might indicate that the regions 
of formation of these lines are not strictly in co-rotation. 
However, the peaks in the periodograms are fairly broad, 
thus these periods are all consistent with a single period- 
icity, between 2.5 and 3 days. 

In the periodogram of Ha, and less conspicuously in 
the periodograms of H/3 and Nai D, there is a broad peak 
at about 5 days that we believe to be related to the flux 
enhancement observed in Ha and H/3 at roughly the mid- 
dle of the time series (Sect. 3.2). We have computed the 
periodograms of the spectral lines, removing the spectra 
affected by this flux enhancement. We confirm our period 
determinations and, in fact, the peak at 5 days vanishes. 

In Fig. |l| we show the average and variance profiles 
of Ha, H/3, Nai D and Hei D3. We also plot a rotation- 
ally broadened photospheric template for the two Balmcr 
lines to indicate clearly the circumstellar contribution 
(Fig. [j]a,b). In both Balmer lines and in Nai D transient 
outflows are present. These appear in the variance profiles 
of these lines as the peaks between —150 and —50 km s — 1 . 
Small mismatches in the wavelength calibration between 
the different sites and the rebinning of the data cause the 
narrow features in the Na I D variance profile at the posi- 
tion of the interstellar absorption lines. 

The average profile of He I D3 is very asymmetric and 
appears to be formed by two components, one at rest ve- 
locity and a redshifted component (Fig. [yd) . The rest com- 
ponent is obtained, in a crude decomposition, by reflecting 
the blue wing of the average profile about zero velocity. 
The redshifted component is isolated by subtracting this 
reflection from the average profile. This component is cen- 
tered at approximately 83 km s _1 . The two components 
are broader than 2 x v sin i and the redshifted component 
extends at least until 200 km s _1 . 



3. Multi-Line Cross-Correlation Analysis 

The spectrum of SU Aur is rich in spectral lines thought to 
originate from the circumstellar environment of the star. 
The Balmer lines, Ha and H/3, the Nai D and He I D3 
lines, due to their different physical properties, probe dif- 
ferent parts of the circumstellar material. Furthermore, in 
these lines mingle several different components: Giampapa 
et al. ( 1993| ) and Johns & Basri ( 19951 ) proposed that the 
profiles of, respectively, Ha and H/3 could be decomposed 
in the contributions of an enhanced chromosphere, and a 
slow and a fast wind components. 

The total time span of our data set allows us to sam- 
ple variability related to rotational modulation. With our 
fine time sampling, we also expect to characterize vari- 
ability on much shorter time scales, for instance arising 
from the inner magnetosphere. We intend to analyse how 
variability propagates through the different spectral lines. 
In our analysis of these complex profiles, we chose the 
cross-correlation function that expresses how the correla- 
tion coefficient between two time series varies with the 
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Fig. 2. Cross-correlation function intensity maps: (a) ACF(R/3), (b) CCF(Ha,H/3), (c) CCF(Nai,H/3) and (d) 
CCF(Hei,H/5)- Each row represents intensity maps for one pair of lines, for 3 time lags: 0, ~ P/2 and ~ P. The pe- 
riod used is approximately 2.75 days. Bright represents high correlation coefficients and dark strong anti-correlations. 
The intensity (grey) scale is not the same in all maps. In each map, the time lag (At) is indicated as well as the 
maximum/minimum correlation coefficients found. 



time delay between them. 

Cross-correlation techniques have been widely used in 
variability studies of active galactic nuclei (White & Pe- 
terson 1994). Smith et al. ( 1999J ) applied this method to 
the equivalent width of the Balmer and Can H & K lines 
for several T Tauri stars, including SU Aur, but their low- 
resolution spectra did not allow to disentangle the differ- 



ent components that mingle in those spectral lines. This 
technique was also previously used, for instance, in Johns 
& Basri ( |l995a ) for the analysis of Ha in several CTTS 
(also SU Aur) and Johns-Krull & Basri ( |I997| ) for the 
analysis of several lines in another TTS DF Tau. 

We calculated the cross-correlation function (CCF) of 
pairs of spectral lines as a function of the time lag r as well 
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as the auto-correlation functions (ACF). In Appendix A 
we define the CCF and describe its implementation for 
our type of time series, as well as the error analysis. 

3.1. Application of CCF Analysis to SU Aur Data 

Using the technique described in Appendix A, we com- 
puted for each pair of spectral lines a cross-correlation 
function intensity map that allows a first glimpse at the 
time lag information contained in the time series. These 
maps were calculated for each pair of spectral lines for lags 
multiple of 0.25 days (6 h). In Fig. g we present the most 
striking examples of the results we obtained. Each column 
in this figure represents the CCF computed at three rep- 
resentative time lags: 0, ~ P/2 and ~ P, where P is the 
rotational period, taken to be approximately 2.75 days. It 
can be seen that the periodicity in the spectral lines is re- 
flected in the behaviour of the CCFs with time lag. The 
CCFs vary from high positive correlation coefficient to 
high negative correlation coefficient in a periodic way (e.g. 
Fig. ga). These maps illustrate the time relation between 
the different spectral lines (Oliveira et al. 2000a). The red 
wing of H/3 seems to vary approximately in phase with the 
red wings of Ha and Nai D (Fig. |b,c respectively), and 
with the whole He I D3 profile (Fig. 2p) . Furthermore, the 
blue and red wings of H/3 seem out of phase (Fig. ||a). 

These kind of maps cannot be used to determine time 
lags with enough precision. Their computation is too time 
consuming to be a viable tool for a finer analysis with bet- 
ter time sampling. Therefore we have decided to choose 
and isolate velocity bins in the spectral lines that are rep- 
resentative of the profiles' variations. 



3.2. Velocity Bins for the CCF 

We selected coherent velocity domains by analysing the 
ACF for each spectral line. This is very well illus- 
trated in Fig. ga for H/3. In this ACF, we can see 
that the H/3 profile is composed of three distinct do- 
mains: [-300:-125] km s -1 , [-125:25] km s" 1 and 
[25:250] km s . Within these velocity domains the profile 
behaves in a concerted way. In each such interval we se- 
lected a representative velocity bin for more detailed CCF 
analysis: —150 km s _1 (the far blue wing), —50 km s _1 
(the slightly blueshifted absorption component) and 
100 km s _1 (the redshifted absorption component). For 
Ha three domains are also defined: [—400: — 125] km s _1 , 
[—125:25] km s _1 and [25:250] km s -1 ; the appropriate 
bins are then —200 km s _1 , —75 km s _1 and 100 km s _1 , 
corresponding to the same regions in H/3. For Nai D the 
blue and red wings seem to vary independently; we have 
thus chosen two velocity bins: —30 km s _1 and 75 km s -1 . 
The profile of Hei D3 varies as a whole within the in- 
terval [—75:150] km s -1 . Still, the Hei D3 profile seems 
composed of two components (Sect. 2.2). Thus, we will 
investigate their behaviour separately, selecting bins at 
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Fig. 3. Intensity variations of the chosen velocity bins. 
From top to bottom: the variations of Ha, H/3, Nai Dl and 
Hei D3, as indicated on the upper corner of each graph. 
Also indicated is the corresponding velocity, on the left 
side. The flux enhancement is clearly visible in the red 
and far blue wings of both Balmer lines. In the slightly 
blueshifted bin of both Balmer lines, from t ~ 12 onwards, 
the effect of the second transient feature can be seen (Sect. 



5). 



and 83 km s _1 . The selected bins for these spectral lines 
are shown in Fig. [j] In Fig. || we plot the intensity vari- 
ations of the velocity bins that we have selected for each 
line: from top to bottom, Ha, H/3, Nai Dl and Hei D3. 

The Ha intensity variations clearly show the flux en- 
hancement seen both in the blue and red wings of the 
profile. The enhancement lasted for a minimum of 2 days, 
and the peak intensity occurs at about t = 9 days. The in- 
crease was ~ 65 % in the blue wing and ~ 72 % in the red 
wing. The central absorption feature (—75 km s _1 ) shows, 
from t ~ 12 onwards, the presence of the second transient 
features (see Sect. 5). In the Ha flux variations there are 
hardly any visible traces of the ~ 3 day periodicity. 

The analysis of the flux enhancement in H3 is more 
complicated as there are fewer data points. Still, the 
flux enhancement can also be seen, causing variations of 
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Fig. 4. Intensity variations folded in time. From top to 
bottom: H/3 at 100 km s _1 , Nai Dl at 75 km s -1 , Hci D3 
at 83 and km s _1 . The intensity variations are plotted 
against: • time, o (time — P), * (time — 2xP). The period 
used was P = 2.63 days. Obviously, the variability of H/3 
has other non-periodic components. Folded in this way, 
the flux variations in He I D and Na I D clearly show their 
periodicity. The lagged behaviour of the variations of the 
three spectral lines is evident. 



~ 22 % in the far blue wing and ~ 37 % in the red wing. 
In the red wing of H/3, the ~ 3 day periodicity can be seen 
and in the slightly blueshifted absorption feature there are 
some indications of this too. Furthermore, the quasi anti- 
phase behaviour of these two parts of the H/3 profile can 
be observed. The variations of the far blue wing of H/3, 
with the exception of the flux enhancement, show no re- 
lation with the variations in the rest of the profile. 

Both components of He I D3 and the red wing of 
Nai Dl, show clear 2.63 days periodicity. The variations 
are not sinusoidal but the periodic deepening of the ab- 
sorption features is clear. The behaviour of the blue and 
red wings of Nai Dl seems anti-correlated. 

In order to clarify and compare more easily the vari- 
ability of the four velocity bins where signatures of period- 
icity were seen, we folded in time the intensity variations 
of H/3 at 100 km s"\ Nai Dl at 75 km s" 1 and Hei D3 
at and 83 km s _1 . This is done in Fig. by plotting 
the intensity variations as a function of time, (time — P) 
and (time — 2xP), starting at to = 6.4 days. The pe- 
riod used was the one determined for Hei D and Nai D: 



P = 2.63 days. We have not done this for any bin in Ha 
because large sporadic changes seem to dominate. In H/3 
there is also some non-periodic variability. The behaviour 
of the other two lines is mainly controlled by the phe- 
nomenon with periodicity of ~ 2.63 days. Very interesting 
in Fig. is that it makes clearly visible the time-lagged 
behaviour between the variability of Hci D3 at and 
83 km s"\ Nai Dl at 75 km s" 1 and H/3 at 100 km s _1 . 
In the remainder of this section we calculate these time 
lags and relate them to accretion and wind flows. 

3.3. Time Lags between Multi-Line Variations 

We compute the CCF for each pair of chosen velocity 
bins with a sampling At = 0.025 day (approximately 
half of the maximum exposure time) for time lags rang- 
ing from —4 to 4 days (for larger time lags the number of 
points in the calculation would be too small) . For each lag 
r the CCF is computed with a different number of pairs 
of points, so the statistical significance of a given value of 
the correlation coefficient varies with the lag. In our anal- 
ysis of the CCFs we will consider only peaks significant 
at more than the 99 % level. We have computed the false 



alarm probability following Bevington & Robinson (1992). 

In Fig. M we show some of the CCFs for the pairs of 
velocity bins that give the most interesting results. In each 
plot, the CCF is represented with error bars that are ob- 
tained from the standard deviation of the Monte Carlo 
simulations. The periodicity of the CCFs is evident and 
is a consequence of the periodic variations in the spec- 
tral lines analysed. Taking for example Fig. 0f, there is a 
central peak of the CCF at about r = 0, and two more 
peaks respectively at ± P. For velocity bins where only 
a weak peak was detected in the period analysis, e.g. the 
blue wing of H/3 and the red wing of Ha, that is reflected 
in the amplitude (i.e. correlation coefficient) and shape of 
their CCFs with respect to the red wing of H/3 (Fig. |]a,c). 
On the contrary, the CCFs between clearly periodic time 
series show sharp peaks with high correlation coefficient, 
as it is seen for the CCF between the redshifted compo- 
nent of Hei D3 and the red wing of H/3 (Fig. ge). 

The shape of the CCFs could be interpreted as mean- 
ing that there are features which appear around a limb of 
the star, last coherently for half a rotation and then dis- 
appear around the other limb. The fact that the "peaks" 
and "valleys" have similar amplitude could indicate that 
features on both sides of the star behave similarly. The ex- 
ception is the case of the two Hei D3 components (Fig. Hf). 

In Table we list the relevant time lags found for each 
pair of velocity bins (only the peak closer to zero was anal- 
ysed in detail). Pairs of velocity bins that did not yield sig- 
nificant lag information or were redundant were discarded. 
We give both the peak and amplitude of the CCF, the me- 
dian of the cross-correlation peak distribution (CCPD) 
and the intervals r ± At67, obtained by interpolation of 
the CCPD (Appendix A). Also listed is the peak's signif- 
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Fig. 5. Cross-correlation as a function of time lag. In each graph the CCF of a pair of selected velocity bins is plotted 
with the error bars obtained from the Monte Carlo simulations; the spectral lines and the velocity bins used are 
indicated. The CCFs were calculated with a sampling of At — 0.025 day. The recurrence of the peaks is a clear 
result of the periodicity present in the spectral lines. The amplitude and shape of the CCF reflect how strong was the 
periodicity detected in the velocity bins. 



icance by means of the false alarm probability (FAP). 

In this table there are positive and negative lags. This 
is a matter of convention in the way the CCF is computed. 
In our case if CCF(l\ ,h) yields a lag r > this means that 
I2 is t in advance of b; conversely if CCF(1 3 ,1 2 ) yields a 
lag r < it means that 1 3 is \t\ in advance of 1 2 . This is 



relevant when several (velocity bins of) different lines are 
compared and sorted in time. 

The far blue and red wings of Ha (respectively —200 
and 100 km s _1 ) correlate well, with a single broad peak 
at approximately zero lag. Also well correlated at zero 
lag with these Ha regions is the far blue wing of H/3 
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Table 1. Time lags found for several pairs of velocity bins. Pairs that do not yield significant lags or provide redundant 
information are not listed. CCF(\i).2) is the cross-correlation function of U and I2. In the first four columns, the spectral 
lines and velocity bins are given. The peak and amplitude of the CCF are listed in columns 5 and 6, the median (r) 
and r ± Atq^ of the cross-correlation probability distribution in columns 7, 8 and 9. The statistical significance of the 
peak by means of the false alarm probability (FAP) is given in the last column. The first six CCFs are listed in the 
order as they appear in Fig. g. 



T50 km s 1 ). This is probably due to the contribution 



of a turbulent broadening mechanism (e.g. Basri 1990; Ed- 
wards et al. |1994| ; Alencar & Basri |2000 ) . The peak of the 
CCF between the red wings of the two Balmer lines is at 
approximately 0.3 days (Fig. ||c), indicating a delay be- 
tween the variability effect on the two lines. The CCF of 
the blueshifted absorption feature in Ha (—75 km s _1 ), 
and the equivalent feature in H/3 (—50 km s _1 ), shows a 
very broad peak at approximately lag r = 0. As a whole, 
the Ha profile correlates less well than H/3 with the other 
lines. This is because the different effects mingle more in 
Ha as it originates in a larger volume. The blue and red 
absorption features in H/3 (—50 and 100 km s _1 respec- 
tively) are anti-correlated at lag r = 0, and show two peaks 
at lag t = ± P/2 (Fig. |a). 

The red absorption features of Nal D correlate well 
with the red absorption feature in H/3 (Fig. ||d) and less 
well with the red wing of Ha. The blue absorption wings 
of Nal D are well correlated with the slightly blueshifted 
absorption features in the Balmer lines. At approximately 
zero lag, the blueshifted regions of these spectral lines are 
anti-correlated with the red wings of the same lines. 

There seems to be a small lag between the variabil- 
ity that affects the two He 1 D3 components, respectively 
and 80 km s _1 (Fig. m). In agreement with this, the 
CCFs of the redshifted absorption features in Na 1 D and 
the two He 1 D3 components also show slightly different 
peaks. However, the other spectral lines, when correlated 
with the He 1 D3 components, yield peaks of the CCF 



that are consistent with a single time lag. Therefore, in 
Table [l] we have only listed the CCF information on one 
of the He 1 D3 components, the component at 80 km s _1 . 
We believe that the time lag detected between the two 
He 1 D3 components is real but small, thus not always 
detectable. Both Hei D3 components are well correlated 
with the red wing of H/3 (Fig. Be) and Nal D. The slightly 
blueshifted absorption features in Ha, H/3 and Nal D are 
anti-correlated at lag r = with He 1 D3 and their CCF 
shows then two peaks at lag r = ± P/2 (Table [l]). 

Between the redshifted absorption component in 
Nal D and the rest absorption component of Hei D3 a 
time lag of —0.075 days was found, indicating the advance 
of the variability in this Hei D3 component. On the other 
hand, the variability in the redshifted Nal D component 
seems to anticipate the variability in the redshifted ab- 
sorption component in Hei D3 by 0.050 days. This time 
lag is at the limit of our time resolution, so we believe the 
variability in these components to be practically simulta- 
neous. The variability of the redshifted absorption com- 
ponent of H/3 is delayed by about 0.28 days with respect 
to the variability in the redshifted Hei D3 component. 

In summary, the absorption features in the red wings of 
H/3 and Na 1 D and the He 1 D3 line correlate strongly over 
a large velocity range. The correlation with the red wing 
of Ha is less clear. Anti-correlated with the red wings of 
H/3 and Nal D, and with Hei D3, are the absorption fea- 
tures in the near blue wings of Ha, H/3 and Nal D. The 
far blue (emission) wings of both Balmer lines are well 
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correlated with each other. We found increasing time lags 
between the rest absorption component of Hci D3, the 
rcdshifted absorption components of Nal D and Hei D3 
and the redshifted H/3 absorption component. 

4. Interpretation: Accretion, Winds and Disk 

4--1- Multi-line Time Lags in SU Aur 

Normally, in stars of spectral type G2 there is no He I D3 
line. The conditions for the formation of the Hei D3 line 
(high temperature or density) are consistent with this line 
being formed at the footpoints of the accretion column 
or slightly above the photosphere, where the kinetic en- 
ergy is released, or by XUV back-warming of the surface 
by an active corona. For most CTTS the He I D3 line ap- 
pears in emission but in SU Aur this line is in absorption. 
We decompose the He I D3 average profile in two compo- 
nents (Fig. |l]d). The component centered at rest velocity 
probably originates in post-shock gas at the base of the 
accretion column, and the redshifted component in the 



accretion column (Najita et al. 200C). We have some indi- 
cation that there is a time-lagged behaviour between the 
two components, and this would also be consistent with 
the two lines being formed in slightly displaced regions. 
But, the fact that we observe the same period (2.63 days) 
and that they vary almost in synchrony confirms that the 
formation of both components is closely linked. 

The absorption component in the red wings of Nal D 
and the He I D3 line vary in an almost synchronized way, 
implying that the Nal D lines are also formed close to 
the star. There is a small delay in the variability observed 
in He I D3 and the variability observed in the absorption 
component in the red wi ng of H/3. Magnetospheric mod- 
els (e.g. Hartmann et al. 1994 ) strongly suggest that the 
redshifted absorption features in the Balmer and metallic 
lines in CTTS are formed in the accretion funnels. We con- 
clude then that the Nal D and H/3 redshifted absorption 
features appear to arise in the inner accretion flow, with 
H/3 forming furthest away from the star. In Ha the situ- 
ation is more ambiguous as there is a large bulk emission 
whose origin is uncertain. Still, the redshifted absorption 
component in Ha is probably also related with the accre- 
tion flow, but at larger distances from the star, and the 
time lags we observe are consistent with this. 

The time lags we determined suggest that a perturba- 
tion is seen in succession in the central absorption com- 
ponent in Hei D3, in the red wing of Nal D and the red- 
shifted component of He I D3, the red wing of H/3 and fi- 
nally the red wing of Ha. We should emphasize two things. 
First, these absorption features do not appear at the high 
velocities characteristic of free-fall. This is consistent with 
viewing SU Aur close to edge-on, i.e. under a small angle 
with the plane of the disk: hence the line-of-sight always 
makes a considerable angle with the base of the accretion 
flow probably at high latitudes on the star. Second, the 
time-lagged behaviour that we observe cannot be due to 



time-variable (non-steady) accretion, because in that case 
the lines that form closer to the disk would be affected 
first and that is the reverse of what we observe. In the 
oblique magnetospheric model, the column density of the 
accreting matter seen in absorption against the star varies 
with the viewing phase of the system: at a certain phase 
we look through the densest parts of the accretion flow, 
whilst half a rotation later the densest parts are no longer 



seen in front of the star (see Johns & Basri 1995b). In 



this scenario no time-lagged behaviour will be observed 
if the magnetosphere co-rotates with the star as a solid 
axi-symmetric body, due to the radial alignment of the 
regions of formati on of the different spectral lines. 

Smith et al. ( 1999 ) found a similar sequence of in- 
creasing time lags for decreasing Balmer series, on one 
night of observations, that they attribute to an accretion 
stream and shock seen as they moved over the limb of the 
star. However, any such occultation effects will be time- 
symmetric: if the ingress of a region of enhanced column 
density is seen subsequentially in spectral lines of different 
excitation conditions, then the reverse sequence is seen for 
the egress of that region. As a result, although the dura- 
tion of eclipse may differ from line to line, the time of 
mid-eclipse will be equal for all spectral lines. Similarly, 
the inner regions of the densest part of the accretion flow 
appear longer in front of the star than the outer regions, 
but the times of mid-transit are the same. 

We propose that the time-lagged behaviour that we ob- 
serve may occur if the anchors of the magnetic field lines in 
the star have a higher rotation rate than those in the inner 
disk. Then the field lines and hence the accretion flow are 
no longer radial and the time symmetry for ingress/egress 
of transit (or occultation) events is broken: mid-transit 
(or mid-eclipse) will occur at different times for spectral 
lines formed at different distances to the star, causing the 
observed time lags. This seems also to be supported by a 
possible increase of modulation period with distance from 
the star (see Sect. 2.2), suggesting that the inner parts of 
the accretion disk are rotating slower than the stellar sur- 
face. Thus there is an external torque working on the mag- 
netosphere that would wind up the magnetic field lines. 
Apparently an equilibrium is established between the ex- 
ternal torque and the magnetic pressure contained in the 
azimuthally distorted magnetic field (see below) . It would 
also hint at a somewhat smaller stellar rotation period of 
P ~ 2.63 days, more easily reconciled with a smaller stel- 
lar radius than the usually adopted R* ~ 3.6 Rq, which 
is also suggested by recent modelling of the spectral en- 
ergy distribution (Oliveira et al. 2000b). 

A time lag r ~ 1.67 day separates the variations in 
the absorption features in the red and blue wing of H/3. A 
slightly bigger time lag (1.77 days) is observed between the 
blue absorption feature in H/3 and the Hei D3 line. Both 
the blueshifted absorption features in Ha and Na I D show 
similar lagged behaviour (time lags of respectively 1.90 
and 1.88 day) with respect to Hei D3. In the scenario of 
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the magnetospheric models, these blueshifted absorption 
features are wind indicators. Thus, the signatures of ac- 
cretion and disk winds are out of phase. Our CCF calcu- 
lations indicate a phase shift of about 240°, but, as it can 
see from Table p} the peaks are broad and their amplitude 
relatively small, so this value is very uncertain. 



4-2. Comparison with Predictions from Magnetospheric 
Models 

We have investigated several magnetospheric models that 
appear in the literature and address the issue of differ- 
ential rotation in magnetosphere and disk. Some mod- 
els that approach the problem of magnetospheric accre- 
tion from the point of view of spin-up/spin-down of the 
young sta rs an d jet formation are reviewed in Goodson & 
Winglee ( 1999J and references therein) . They all consider 
the case of strong magnetic fields (i.e. there is disk trun- 
cation) and the physical distinction between them comes 
from the effective value of magnetic diffusion in the disk. 
In these models, the different rotation rates of the star 
and the disk cause the field lines that permeate the disk 
to become azimuthally twisted and a toroidal component 
of the magnetic field is created. The ability of the disk 
to counteract this toroidal component is described by the 
effective magnetic diffusivity. 

For the most diffusive flows, no wrapping up of the 
magnetic field lines occurs because the field lines can con- 
tinuously slip through the disk. This is the case for the 



Ghosh & Lamb fll978| ) and the Shu et al. Ql994| ) models. 
According to Ghosh & Lamb all the field lines are closed 
and they thread the disk up to distant radii. In the Shu et 
al. model the closed field lines thread only a small, nearly 
uniformly rotating, region of the disk and outside this re- 
gion there are open magnetic field lines. 

For low magnetic diffusivity, differential rotation be- 
tween the disk and the star makes the magnetic field lines 
wrap up and open via a mechanism known as inflation 
(Lovelace et al. 1995 ). The inner edge of the disk under- 
goes oscillations and magnetic reconnection occurs caus- 
ing accretion episodes ; no s teady configuration is possible 
(Goodson & Winglee |1999| ). 

For intermediate diffusivity, once inflation opens the 
outer magnetic field lines, a steady state can be attained. 
The inner magnetic field lines remain closed and approx- 
imately in co-rotation with the star. The plasma in the 
disk can move radially inwards until it reaches the point 



where steady accretion occurs (Lovelace et al. 1995). In 
this last scenario, it is plausible that the closed magnetic 
field lines are azimuthally deformed. 

The time-lagged behaviour that we observe in SU Aur 
suggests a quasi-steady configuration, in which the mag- 
netic field lines are azimuthally curved as a result of the 
magnetic field line anchors being dragged by, but slipping 
through the inner disk that is rotating somewhat slower 
than the magnetic field line anchors on the stellar surface. 



The misalignment of rotation and magnetic axes intro- 
duces, at certain phases, an extra magnetic pressure com- 
ponent in the plane of the disk, that eases the loading of 
the closed magnetic field lines with matter from outside 
the co-rotation radius. 

But how do the outflows relate with this proposed sce- 
nario? Our proposed scenario includes open field lines that 
can drive wind flows (as suggested by Shu et al. 1994). We 



have found that the signatures of accretion and disk winds 
are about 240° out of phase. This agrees with the sugges- 
tion of Johns & Basri ( 1995b| ) (confirmed by Petrov et al. 



1996 ) if their picture of the oblique rotator is refined to 



include azimuthally curved magnetic field lines. 



5. Transient Blueshifted Absorption Features 



5.1. Velocity Variations 



Two transient absorption features were detected in the 
blue wings of the profiles of Ha, H/3 and Nai D (Un- 
ruh et al. 1998 ). We did not pursue a multi-comp onent 
analysis of these lines (as done in Giampapa et al. 1993, 



Johns & Basri 1995b or Oliveira & Foing 1999) due to 
the optical thickness of the emission and the complexity 
and dynamics of the profiles. We could however isolate the 
transient features in the spectra by applying an adequate 
filtering procedure. The velocity position of the features is 
determined by detecting the minimum (in a chosen veloc- 
ity interval) in the difference spectra with respect to the 
median-filtered spectra. In Fig. o, the velocity evolution 
of these features is plotted (features A and B). There is 
no considerable change in the width/depth of the features, 
but these measurements are very difficult to carry out due 
to the intrinsic variability of the underlying profile. 

Feature A persisted for about 1.6 days. It drifted from 
approximately —140 to —90 km s _1 in Ha and from —110 
to —70 km s _1 in Nai Dl. In H/3 we cannot identify the 
velocity position of the feature unambiguously and so the 
minimum detection was more uncertain. However, there 
are indications that |v(Ha)| > |v(H/3)| > |v(Nai)|. The 
velocity evolution in the three lines follows parallel tracks 
in time. The onset of the feature A is simultaneous with 
the flux enhancement observed in Ha and H/3 (Sect. 3.2). 

Approximately 4 days after the onset of feature A, an- 
other transient feature was detected in Ha and H/3 but 
no counterpart was seen in the Nai D lines (feature B in 
Fig. 0) . Feature B may therefore originate in regions higher 
in the magnetosphere that contribute little to the Nai D 
line formation. In Ha this feature drifted from —140 to 
—80 km s _1 and seems to be decaying after 3 days when 
our time series ends. In contrast to what we observe for 
feature A, the velocity position of feature B in H/3 is prac- 
tically the same as in Ha. We stress that feature B is more 
conspicuous in H/3. In the case of feature B it is rather 
difficult to be certain if it is an absorption or emission 
transient, or a combination of both. 
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Fig. 6. Velocity evolution of two drifting absorption fea- 
tures in the line profiles of • Ha, o H/3 and x Nai Dl. For 
Nal D2 the velocity position of the feature is similar to 
Nai Dl. Feature B was not detected in Nai D. Feature A 
lasts for at least 1.6 days and feature B for 3 days. The 
velocity position of the features is determined by detect- 
ing the local minimum in the line profiles. Simultaneously 
with feature A, an intensity enhancement was observed in 
Ha and H/3, affecting the whole profiles. 



5.2. Interpretation: Sporadic Mass Ejections 

We propose that the first transient absorption feature 
might be plasma ejected upward, that condenses and after- 
wards decelerates under the influence of gravity, causing 
the observed velocity evolution. We have investigated sev- 
eral scenarios that could cause such velocity evolution. A 
blob of matter moving at constant velocity, in which case 
the velocity variation would be just the result of angu- 
lar projection, is ruled out because it could not cause the 
velocity difference between the lines. The same problem 
is faced when postulating a blob that departs from the 
stellar surface with the escape velocity and then deceler- 
ates. The observed offset between the initial velocities of 
the lines indicates radially elongated material that under- 
went an acceleration outwards by an as yet unidentified 



mechanism, before it started to decelerate. For the second 
transient feature the velocity evolution shows no evidence 
of any initial acceleration. Still, as for the first feature, 
the velocity evolution strongly suggests deceleration even 
though the dispersion is very large. 

We investigated several types of solar-like phenomena 
that could have similar characteristics. Post-flare loops or 
surges are out of consideration because, although the mat- 
ter is initially accelerated and condenses, it then falls back 
to the solar surface and would thus appear as redshiftcd 
absorption features, not blueshifted as we observe. Such 
redshifted absorption components have been d etecte d on 
other stars though, e.g. H R 109 9 (Foing et al. |l994[ ) and 
AD Leo (Houdebine et al. 1993 ) and were proposed to be 
the analogue of post-flare loops. 

If the deceleration measured for feature A is inter- 
preted as gravitational deceleration of material, then a fit 
to that deceleration rate indicates a location of the decel- 
erating matter at approximately 7 to 8 R*. It is difficult to 
explain a column of material progressing at such distances 
from the star, specially without knowledge of its previous 
evolution. 

Feature A might be the line-of-sight absorption of ma- 
terial of a large scale coronal mass ejection, in a config- 
uration sometimes called a halo mass ejection. Similar to 
what is believed to happen in the Sun and flare stars (Fo- 
ing 1989] ; Houdebine et al. 1993), a disruption of the mag- 
netic field lines, perhaps near the truncation radius, might 
cause a "forward" mass ejection, away from the star, that 
cools and condenses and is observed as transient absorp- 
tion features, and a "backwards" release of matter/energy, 
that can hit the star and cause the Balmer brightening 
we observe. In the case of SU Aur, these ejections clearly 
point at a magnetospheric connection between the star 
and the disk where magnetic stresses and energy build 
up and are released sporadically. To reproduce the time 
evolution of the flare, as seen in the Balmer brightening, 
(peak at day 9.1, Fig. ||) and of the decelerated absorp- 
tion component (start at day 8, Fig. ||), we suggest that 
the primary instability and energy release did not occur at 
the stellar surface or in the lower atmosphere, but rather 
at high altitude in the interacting magnetosphere between 
the star and the disk. For the second absorption feature, 
a similar scenario would apply, but without a visible de- 
tection of a flare. 

The events that we observed are very likely triggered 
by magnetic field disruption, perhaps in the context of 
sheared magnetic field lines as described in the previous 
section. The energy released during these events, in combi- 
nation with more continuous energy release from magnetic 
field reconfiguration, might explain the especially high X- 
ray luminosity of SU Aur (Skinner & Walter 1998). 
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6. Summary 

The circumstellar environments of CTTS are extremely 
interesting for studying the interaction between the star 
and the disk. The CTTS line profiles are very complex 
and diverse, as both accretion and outflows play a role 
in their formation. In this paper we try to disentangle to 
some extent these effects, by analysing Ha, H/3, Nai D 
and He I D3, in an attempt to constrain the geometry and 
dynamics of the magnetosphere. 

The He I D3 line in SU Aur has two components. Their 
variability is slightly time-lagged, indicating spatially dis- 
placed regions of formation. The component at rest ve- 
locity can be formed close to or at the stellar surface, for 
instance in the post-shock region, and the redshifted com- 
ponent in the accretion column. We propose that the red- 
shifted component of He I D3, the red wings of Nai D and 
the redshifted absorption features in H/3 and Ha are all 
formed in the accretion funnel, at different distances from 
the star according to their opacities and excitation condi- 
tions. We find different time lags between these lines, with 
the variability of the lines closer to the star preceding the 
variability in the outer spectral lines. We argue that such 
behaviour cannot be caused by non-steady accretion, nor 
by simple occultation effects in an oblique magnetospheric 
model. We propose that such time-lagged behaviour can 
only occur if the magnetic field lines in the oblique mag- 
netosphere have an azimuthal component. 

We find that the near blue wings of Ha, H/3 and Nai D 
(wind indicators) are out of phase with the accretion in- 
dicators, though not exactly in anti-phase. 

We observe two transient absorption features in the 
blue wings of Ha, H/3 and Nai D. The velocity evolution 
of the first feature indicates matter to be initially accel- 
erated (different velocities for the different lines). Subse- 
quent evolution suggests a deceleration or a projection 
effect. We propose that these features are the signature of 
plasma outflow events, powered by the disruption of the 
distorted magnetic field lines. 
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Appendix A: Cross-Correlation Function Analysis 



A.l. Definition and Computational Method 

Let us assume two stationary time series Li(t) and L2(t) 
discretely sampled N times at times t,, with At = t»+i — t, 
for all values 1 < i < N— 1. The cross- correlation function 
(CCF), computed at intervals r (or time lags) which are 
integer multiples of the sampling interval At, is defined as 
follows: 

CCF(r) = i £ [L 1 fa)-L 1 ][L 2 (t t -r)-L 2 ] 



cr 1 a 2 



where Lk and Ck are respectively the mean and standard 
deviation of each time series Lk(tj). The obvious problem 
is that our time series is unevenly spaced, so this defini- 
tion has to be adapted. 

There are two methods to circumvent this difficulty: 
the interpolation method and the discrete correlation func- 



tion. In the first method (Gaskell & Peterson 1987), the 



CCF is calculated twice for each time lag r: the observed 
points Li(tj) are correlated with the interpolated values 
of L 2 (tj — t) and the observed points L 2 (tj) are correlated 
with the interpolated values Li(tj+r). The two results are 
then averaged, assuming there is no reason to prefer one 
interpolation over the other. The second method, the dis- 
crete correlati on fun ction (DCF) was introduced by Edel- 
son & Krolik ( 1985 ). In this method the cross-correlation 
function is computed at each time lag using only real data 
points separated by r. The resulting correlation function 
is then binned on intervals St so that the calculated value 
at time lag r is an average over the interval r ± St/2. 

These methods were found to be in good agreement for 
well sampled data sets with a large number of data points 



(White & Peterson 1994, henceforth WP, and references 
therein). The DCF method is more general in the sense 
that it makes no assumptions on the nature of the time 
series that we are sampling. When using the interpolation 
method, it is assumed implicitly that the sampled time 
series vary smoothly on the time scales that correspond 
to the intervals between observations. Still, the DCF will 
require more points than the interpolation method to give 
a meaningful result and, in the limit of poor sampling, 
will more easily fail to give a significant lag. Therefore, 
for our data set we decided that the interpolation method 
is more adequate and hence used this in the calculations 
presented in this paper. 

A. 2. Optimization Techniques 

When calculating the CCF for several types of test func- 
tions with the same time sampling as our data points, we 
realised that the calculations yielded CCF peak values 
above unity. There was clearly a normalization problem 
in the calculation of the CCF. For stationary time series, 
the mean and variance of the sample are invariant in time; 
only in this situation it is valid to calculate Lk and Ck using 
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the entire time series. We have concluded that this is not 
always a reasonable assumption: in fact, in our data set 
there was a strong flux enhancement in the spectral lines 
(Sect. 3.2), so clearly these statistical quantities change 
in time. Thus, we concluded that it was more appropriate 
to compute Lk and Ck for each lag r, with only the data 
points effectively used for the computation of CCF(t). 
Therefore, the calculation of CCF(t) reduces to the cal- 
culation of the linear correlation coefficient at each lag r. 

We address now the question of how the interpolation 
is actually performed. First, following WP, we do not in- 
terpolate beyond the limits of the time series: we exclude 
from the calculation of CCF(t) points where t» — r < ti 
or t j — t > tN ■ We also use linear interpolation of the data 
points (WP). As it is mentioned above, the interpolation 
assumes that the times series vary smoothly between the 
sampled data points. In our time series the extreme data 
points are separated from the rest of the observations by 
more that 1.5 days. We have computed the CCF both in- 
cluding and excluding these points and we concluded that 
the CCF peaks are much sharper in the last case. So we 
have excluded those data points from CCF calculation, 
even though we still use them for the calculation of the 
interpolated series. We should point out that the number 
of points used is only decreased to N— 4. 

As we have mentioned the sampling of H/3 is worse: 
we have only 72 H/3 spectra against 126 for the remain- 
ing lines. The calculated CCF(t) is an average quantity 
if there is no reason to favour one of the interpolations. 
But if we are interpolating H/3 and Ha, for example, it is 
clear that the interpolation of the Ha series is better be- 
cause it is computed with more data points. Thus, when a 
CCF that involves H/3 is computed we will consider only 
one term, the one where the other time series is interpo- 
lated. In the example above we consider only the terms 
where the interpolated Ha series is cross-correlated with 
the observed H/3 series. 

A. 3. Error Analysis 



Gaskell & Peterson ( |1987|) propose an analytical formula 
to compute the uncertainty of the time lag r. WP and 
Maoz & Netzer ( 1989[ ) found that this formula can severely 
underestimate the uncertainties. Accordingly, we have de- 
cided to use Monte Carlo simulations to estimate these 
errors. We assume that the noise in our time series has 
a gaussian distribution. We use the standard deviation of 
the continuum as a conservative estimate of the noise in 
the spectral lines. We computed the CCF with a sam- 
pling step of At — 0.025 day (approximately half of the 
maximum exposure time in the time series). We used 1000 
trials for the simulations. We then computed the empirical 
cross- correlation peak distribution (CCPD) (Maoz & Net- 
zer 1989), by simply counting how many times the peak 



occurred at each time lag. The peak of the CCPD is the 
most likely time lag and its width is the uncertainty of this 



estimate. Given the asymmetry of the CCPD, we follow 
WP and define the quantity ±At67, for which 2/3 of the 
trials give a lag in the interval r± Atqi with relation to the 
median of the distribution. When the CCPD is a gaussian 
At67 corresponds to the standard deviation. The quantity 
Ar 6 7 is derived directly from the cumulative CCPD. 
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